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I.  INTRODUCTION 


In  this  paper  we  derive  the  governing  equations  for  a  lumped  parameter 
model  of  a  regenerative  liquid  propellant  gun  and  discuss  the  numerical 
Implementation.  Morrison  et.  al.  have  covered  the  background  and  possible 
applications  of  liquid  propellant  guns.^  Here  the  emphasis  Is  on  the 
derivation  of  a  numerical  model. 

2 

For  solid  propellant  guns,  there  Is  a  wide  range  of  models. 

Precomputer  models  Involved  approximations  that  allowed  an  analytic, 
closed-form  solution,  or  else  a  simple  numerical  solution  leading  to 
tables.  The  advent  of  computers  allowed  the  direct  numerical  solution  of 
ordinary  differential  equations  describing  the  Interior  ballistics  on  a  case 
by  case  basis.  More  recently,  one-dlmenslonal  and  even  two-dimensional 
models  have  been  developed.  Involving  the  solution  of  partial  differential 
equations.  All  of  these  procedures  can  give  accurate  results.  The  choice 
of  a  model  depends  on  the  detail  desired  In  the  solution. 

For  regenerative  liquid  propellant  guns,  much  less  work  has  been 
done.  Pagan  and  Izod^  have  developed  an  ordinary  differential  equation 
model.  This  code  assumes  a  particular  type  of  regenerative  gun,  and  has  a 
rather  complicated  model  for  the  liquid  Injection  Into  the  combustion 
chamber.  Cushman^  has  written  a  model  that  has  been  extensively  used  at 
General  Electric.  The  model  normally  used  at  the  BRL  Is  due  to  Paul 


^W.F.  Morrison,  J.D.  Knapton,  and  G.  Klingenberg,  "Liquid  Propellants  for 
Gun  Applications,”  Proceedings  of  the  Seventh  International  Symposium  on 
Ballistics,  The  Hague,  The  Netherlands,  April,  1983. 

^P.G.  Baer,  "Practical  Interior  Ballistic  Analysis  of  Guns,"  Progress  In 
Astronautics  and  Aeronautics,  Interior  Ballistics  of  Guns.  (H.  Krier  and  M. 
Summerfleld,  ed.),  Vol  66,  1979. 

^G.  Pagan  and  D.C.A.  Izod,  "Regenerative  Liquid  Propellant  Gun  Modelling," 
Proceedings  of  the  Seventh  International  Symposium  on  Ballistics,  The  Hague, 
The  Netherlands,  April,  1983. 

^P.G.  Cushman,  "Regenerative  Liquid  Propellant  Gun  Simulation  User's 
Manual,"  GE  Report  84-POD-004,  December,  1983. 
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Gough. ^  This  code  uses  a  lumped  parameter  representation  (ordinary 
differential  equations)  for  the  liquid  reservoir  and  combustion  chamber. 

The  gun  tube  Is  represented  by  partial  differential  equations.  This  code 
has  a  number  of  options,  Including  a  traveling  charge  projectile.  Since 
both  the  design  of  the  regenerative  liquid  gun  and  the  assumptions  made  In 
the  model  are  under  development,  additional  options  are  periodically  added 
to  the  code. 

The  purpose  of  this  report  is  twofold.  First  we  develop  a  complete 
lumped  parameter  code  (only  ordinary  differential  equations)  to  describe  the 
regenerative  liquid  propellant  gun.  This  code  Is  not  as  general  as  that 
developed  by  Gough.  However,  It  runs  much  faster,  which  allows  parametric 
studies  to  be  made  easily.  Also,  the  code  Is  written  as  far  as  possible  in 
a  modular  form.  Additional  options  can  be  added  with  a  minimum  amount  of 
effort.  As  new  gun  designs  are  developed,  or  new  assumptions  about  the 
behavior  In  the  gun  are  tried,  these  can  be  tested  fairly  easily. 


II.  THE  REGENERATIVE  LIQUID  PROPELLANT  GUN 

Figs.  1  and  2  show  regenerative  liquid  propellant  guns  with  an  In-line 
piston.  Another  possible  design  is  the  reverse  annular  gun,  where  the 
piston  Is  wrapped  around  the  barrel  and  moves  In  the  direction  of  the 
muzzle.^  The  reverse  annular  gun  includes  an  intermediate  combustion 
chamber.  While  we  will  only  consider  in-line  pistons,  in  order  to  be 
consistent  with  the  notation  used  by  Gough,  the  regions  of  the  gun  will  be 
denoted  as  regions  1,  3,  and  4. 

Region  1  Is  the  liquid  propellant  reservoir.  The  monopropellant  Is 
pumped  into  the  reservoir  at  the  beginning  of  the  firing  cycle.  A  primer  Is 
ignited  In  the  combustion  chamber  (region  3).  As  the  chamber  is 
pressurized,  the  piston  Is  forced  back.  Because  of  the  piston  area 


^P.S.  Gough,  "A  Model  of  the  Interior  Ballistics  of  Hybrid  Llquld-Propellant 
Guns,"  Final  Report,  Contract  DAAK11-82-C-0154,  PGA-TR-83-4,  September 
1983. 
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Figure  1.  A  Regenerative  Liquid  Propellant  Gun 


Figure  2.  A  Regenerative  Liquid  propellant  Gun  with  an  Annular  Piston 


differential  between  the  two  regions,  the  liquid  pressure  is  higher  than  the 
combustion  chamber  pressure.  So  liquid  propellant  is  forced  through 
circular  holes  in  the  piston  face,  and  ignites  and  burns  in  the  combustion 
chamber.  Eventually,  the  pressure  pushes  the  projectile  down  the  gun  tube 
(region  4). 

In  most  problems  involving  liquids,  the  liquid  can  be  considered 
incompressible.  Because  of  the  high  pressures  in  a  gun,  this  is  not  a  good 
approximation,  and  the  compressibility  of  the  liquid  propellant  will  be 
considered. 

In  the  gun,  it  is  possible  that  gas  from  the  combustion  chamber  can  be 

forced  into  the  liquid  chamber.  This  is  most  likely  to  occur  at  early 

times,  when  the  primer  causes  the  initial  pressure  rise  in  the  combustion 
chamber.  Because  of  the  inertia  of  the  piston,  there  may  be  a  delay  before 
the  liquid  pressure  also  rises.  The  modeling  of  bubbles  in  a  compressible 
liquid  is  complicated,  so  we  will  assume  that  there  is  no  mass  flow  into 

region  1.  In  the  experimental  fixture  grease  is  used  to  block  the  holes  so 

that  flow  into  the  liquid  reservoir  is  less  likely. 

The  orifice  flow  through  the  vents  in  the  piston  is  quite 
complicated.  Here  we  assume  that  the  piston  face  is  infinitely  thin,  and 
ignore  the  behavior  in  the  orifice.  The  mass  of  the  piston  will  be 
considered. 

The  piston  design  considered  here  was  used  in  early  test  fixtures,^  ^ 
but  is  not  convenient  for  actual  guns.  Since  the  vent  holes  must  be  plugged 
before  each  shot,  rapid  fire  is  impossible.  Other  piston  designs  will  be 
considered  below. 


^W.F.  Morrison,  M.J.  Bulman,  P.G.  Baer,  and  C.F.  Banz,  "The  Interior 
Ballistics  of  Regenerative  Liquid  Propellant  Guns,"  1984  JANNAF  Propulsion 
Meeting,  New  Orleans,  LA,  CPIA  Publication  390,  February,  1984. 

^P.G.  Baer  ,  C.F.  Banz,  I.W.  May  and  W.F.  Morrison,  "A  Propulsion  System 
Comparison  Study  for  the  120-MM  Anti-Armor  Cannon,"  1984  JANNAF  Propulsion 
Meeting,  New  Orleans,  LA,  CPIA  Publication  390,  February,  1984. 
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The  behavior  in  the  combustion  chamber  is  very  complicated.  A  primer 
is  ignited  in  the  chamber,  leading  to  a  gradual  pressure  rise.  Liquid  jets 
are  forced  out  of  the  vent  holes.  These  jets  may  break  up  into  droplets 
because  of  hydrodynamic  forces  between  the  gas  and  the  liquid  or  due  to 
impact  on  the  wall  of  the  chamber.  The  droplets  formed  may  break  up  further 
or  coalesce.  The  propellant  will  eventually  ignite,  and  may  burn  as 
individual  droplets  or  as  an  envelope  flame.  These  processes  are  not  well 
understood  at  low  pressure,  and  essentially  no  work  has  been  done  at  gun 
pressures.  Therefore,  simplifying  assumptions  are  necessary. 

As  a  first  approximation,  we  assume  that  the  liquid  burns  instan¬ 
taneously  as  soon  as  it  enters  the  combustion  chamber,  releasing  all  its 
energy.  The  combustion  chamber  is  treated  as  a  homogeneous  region. 

The  primer  is  assumed  to  be  the  same  liquid  as  the  propellant. 

Normally  the  combustion  of  the  primer  is  not  followed  in  detail.  Instead 
the  initial  pressure  in  the  combustion  chamber  is  an  input  parameter.  Then 
the  amount  of  propellant  needed  to  produce  this  pressure  is  calculated,  and 
this  is  taken  as  the  mass  of  the  primer.  The  piston  and  projectile  are  not 
allowed  to  move  during  the  burning  of  the  primer.  Later  a  somewhat  more 
complicated  primer  option  will  be  discussed. 

There  is  some  evidence  that  liquid  accumulates  in  the  combustion 
chamber  at  the  beginning  of  the  firing  cycle,  which  changes  the  behavior  of 
the  pressure  rise.  Later  we  introduce  finite  burning  rates  for  the  liquid 
propellant,  in  an  attempt  to  model,  at  least  crudely,  the  effects  of 
possible  propellant  accumulation. 

Last,  we  have  the  behavior  in  the  gun  tube  itself.  In  the  Gough  code, 
this  region  is  modeled  by  one-dimensional  partial  differential  equations. 
This  is  necessary  because  of  the  traveling  charge  option  included  in  this 
code.  We  will  consider  this  region  as  a  lumped  parameter  region,  using  the 


Lagrange  pressure  dlscrlbutlon.  This  has  been  successfully  used  Co  model 
the  behavior  of  solid  propellant  guns.^ 

Below,  we  derive  the  governing  equations  for  the  simple  lumped 
parameter  model  just  described.  Later,  additional  options  will  be 
considered. 


III.  GOVERNING  EQUATIONS 

With  the  above  assumptions,  the  regenerative  gun  behavior  can  be 
modeled  by  13  ordinary  differential  equations  ([1]  to  [13]). 

The  equations  governing  the  piston  motion  can  be  derived  simply.  Let 

V2  be  the  volume  of  the  fuel  chamber,  and  Aj  be  the  area  of  the  fuel  side  of 

the  piston  (this  area  Includes  the  vent  holes).  Let  be  the  initial 

volume  of  the  liquid  chamber.  Let  S-„  be  the  piston  travel  (to  the  left) 

ps 

and  be  the  piston  velocity.  Then 


^  =  '^lo  -  Sps  ^  • 


Taking  the  derivative. 


■‘1  dr--V*i- 


The  acceleration  of  the  piston  equals  the  force  on  the  piston  divided 
by  the  mass  of  the  piston,  that  is 


®J.  Corner,  Theory  of  the  Interior  Ballistics  of  Guns.  Wiley,  New  York, 
1950. 


■**  h"*  kV  U*  .i'"  I 
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dV  g 

"  Pi 

ps 


where  A^  is  the  area  of  the  combustion  side  of  the  piston,  is  the  area  of 
the  vent  holes,  p^  is  the  pressure  in  the  liquid  chamber,  p^  is  the  pressure 
in  the  combustion  chamber,  and  M  is  the  mass  of  the  piston.  The  quantity 
••  10'  gm/s-cm-MPa  is  a  conversion  constant  to  put  the  acceleration  in  the 
desired  units  of  cm  /s.  By  assumption,  the  pressures  are  constant 
throughout  the  two  chambers.  Then  the  equation  governing  the  piston  travel 
is 


[3] 


V 

ps 


(4) 


As  the  piston  moves  to  the  left,  liquid  propellant  is  forced  through 
the  vent  holes.  This  is  governed  by  Bernoulli's  equation^ 


where  is  the  density  of  the  propellant  in  region  1,  Vj  is  the  velocity 
of  the  fluid  on  the  liquid  side,  and  v^  is  the  velocity  of  the  injected 
fluid  on  the  combustion  chamber  side  of  the  piston.  In  using  this  equation, 
the  flow  in  the  orifices  is  assumed  to  be  frictionless  incompressible 
flow.  The  potential  energy  due  to  height  differences  has  been  left  out  of 
the  equation. 


^W.F.  Hughes  and  J.A.  Brighton,  Fluid  Dynamics,  Schaum's  Outline  Series, 
McGraw-Hill,  New  York,  1967. 
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Since  the  liquid  in  region  1  is  a  homogeneous  mixture,  we  take  Vj  to  be 
zero.  Eq.  (5)  can  be  rewritten  as 


Pj  V3 


2pj  (Pj  -  P3)  . 


The  mass  flow  rate  through  the  orifices,  denoted  by  “13*  equals  P^v^A^  . 
There  is  usually  some  narrowing  or  constriction  of  the  jet.  This  is  taken 
into  account  by  introducing  a  discharge  coefficient  Cjj.  Then 


*13  *  S  \  -  P3^  •  Pi  >  P3  • 


P,  <  P 


1  "  *^3  * 


For  short  circular  orifices  at  atmospheric  pressure,  the  discharge 
coefficient  should  be  about  0.61  to  0.63.^*^  The  orifices  are  plugged  at  the 
beginning  of  the  firing  cycle,  so  that  gas  will  not  be  forced  into  region 
1.  So  the  mass  flux  is  set  to  zero  if  p^  <  p^. 


Now  let  Mj  be  the  mass  of  the  liquid  in  region  1.  By  conservation  of 


dt  “  "  “l3 


The  mass  is  equal  to  the  density  times  the  volume  of  region  1.  Rearranging 
eq.  (8), 


^W.  Kaufman,  Fluid  Mechanics,  McGraw-Hill,  New  York,  1963. 


(9) 


dp 


(4J 


_1 

dt 


Pi  dV^ 


m 


dt 


n 


This  equation  becomes  singular  as  the  volume  approaches  zero.  This  problem 
will  be  discussed  below  in  the  section  on  the  numerical  solution  of  the 
equations. 

The  bulk  modulus  K  of  a  compressible  fluid  is  defined  by^^ 


P 


1 


K  , 


(10) 


where  the  subscript  T  indicates  that  the  derivative  is  evaluated  with  the 
temperature  constant.  The  liquid  is  assumed  to  be  isothermal.  The  bulk 
modulus  does  vary  with  pressure.  For  lack  of  better  information,  we  assume 
a  simple  linear  form 

K  =  Kj  +  K2  Pi  .  (11) 

Combining  eq.  (10)  and  (11),  and  solving  for  Pi , 

K 

Pi  =  [(kPi)  ^  -  KjJ/K^  ,  (12) 

where  k  is  a  constant.  Let  p^  be  the  density  at  p*0.  Usually  the  density 
is  known  at  p  «  1  atm. ,  but  compared  to  gun  pressure  ranges  1  atm.  is 
essentially  zero.  Then  solving  for  k  and  substituting  into  eq.  (12) 


(13) 


^^R.W.  Fox  and  A.T.  McDonald,  Introduction  to  Fluid  Mechanics ,  2nd  ed., 
John  Wiley  and  Sons,  NY,  1978. 


Taking  derivatives, 


l!i  -  ■'i  /!iY2 

'  "i  \y 

and  combining  eq.  (11)  and  (13)  this  can  be  rewritten  as 


f!i  .  iL  ^ 

dt  ~  Pj  dt  * 


(] 


This  can  be  expressed  in  terms  of  the  speed  of  sound  in  the  liquid  c^,  wl 

1  9 

is  given  by^ 


or 


C 


Then  eq.  (15)  becomes 


dp,  c  ^  dp 

151  5^ 


Now  consider  the  combustion  chamber.  By  analogy  with  eq.  [1], 


Band,  Introduction  to  Mathematical  Physics.  D.  Van  Nostrand  Company 
Princeton,  NJ.,  1960. 
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(19) 


ps  "3 


The  equation  for  conservation  of  mass  is 


(20) 


where 


is  the  mass  flux  from  region  3  to  region  4. 


It  follows  that 


,7,  !3  ^  ^3  -  "34 

dt  “  ■  dt  “v 


(21) 


Finally,  we  need  an  equation  for  the  pressure  in  region  3.  This  will 
be  based  on  the  Noble-Abel  equation  of  state® 


P3  -  P3  RT^/Wd  -  p^b)  , 


(22) 


where  T^  is  the  temperature  in  region  3,  R  is  the  universal  gas  constant,  W 
is  the  molecular  weight  of  the  gas,  and  b  is  the  covolume.  This  can  also  be 
trritten  as 

P3  -  P3(y  -  1)  T3/(1  -  P3b)  ,  (23) 


where  c^  is  the  specific  heat  at  constant  volume  and  y  is  the  specific  heat 
ratio  Cp/Cy.  It  will  again  be  useful  to  determine  the  speed  of  sound® 


20 


where  s  is  the  entropy.  Taking  the  partial  derivative  of  eq.  (23),  after 
some  algebra 


®P3/s  ^3^^  "  ^^^3^ 


where  we  have  used  the  fact  that 


3T3\  (y  -  1)13 
^/s  '  • 


Combining  eq.  (24)  and  (25)  results  in 


J-^1 


'3  ■  P3(l  -  P3b)  • 


The  energy  equation  for  a  homogeneous  gas  can  be  written  as^ 


dT3  /31nV3\  3P3 

^3  ^p  ~di  “y3lnT3jP3  IT  *  ^3 


where  the  spatial  terms  and  the  viscosity  have  been  ignored.  The  quantity 

A 

V  “  l/p3*  Using  the  Noble-Abel  equation. 


^^G.N.  Lewis  and  M.  Randall,  Thermodynamics.  McGraw-Hill,  New  York,  1961. 

Bird  and  W.E.  Stewart,  Transport  Phenomena.  John  Wiley  and  Sons,  New 
York,  I960. 


and  eq.  (28)  becomes 


3T.  dp 

'3  ■  "“s’  dT  ■"  ®3 


(30) 


The  source  term  is  due  to  the  energy  release  when  the  liquid 
combusts.  In  general. 


S 


3 


Zh.  R.  M. 

Ill 


(31) 


where  h^  is  the  enthalpy  of  the  quantity  i,  is  the  rate  of  production  of 
the  quantity  i,  and  is  the  mass  of  the  quantity  i.  The  rate  of 
production  of  gas  in  region  3  is  mj^/V^*  since  the  liquid  combusts 
instantaneously  when  it  enters  the  combustion  chamber.  The  outflow  is 
irrelevant  since  it  does  not  cause  a  change  in  enthalpy. 

Again  from  the  Noble- Abel  equation. 


dT 


_3 

dt 


(1  -  bp^) 


dp 


3 

dt 


!3 

^3 


(32) 


Substituting  into  eq.  (28)  and  doing  the  algebra 


18] 


dp^  dp^ 

dt  *  g  dt  ^  V_  -  bM_ 

“o  3  3 


*34  >0 


(33) 


Normally,  gas  will  flow  from  the  combustion  chamber  into  the  gun 
tube.  But  if  the  gas  flow  from  region  3  to  region  4  (^^4)  is  negative, 
there  will  be  an  additional  source  term.  In  this  case 
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1^3  .  '3_  dp, 

r.  iT”KM 


dt  g  df 
®o 


V3  -  bM3 


“34<\  - 

-  b 


The  enthalpy  of  the  liquid  is  given  by 


hj  =  Cl  +  P^/Oj  » 


where  e^^  is  the  chemical  energy  of  the  liquid.  This  is  obtained  from  closed 
bomb  tests.  The  enthalpy  of  the  gas  is  given  by 


^3  ~  ^v  '^3  ^3^  ^3  ~  ^p  '^3  * 


Finally,  we  consider  the  gun  tube  (region  4).  The  volume  has  the 
standard  equation 


dE-  =  '^pj  ^  • 


where  Vpj  is  the  velocity  of  the  projectile  and  is  the  area  of  the  gun 

tube.  To  model  this  as  a  lumped  parameter  region,  we  assume  a  Lagrange 

Q 

pressure  distribution.”  That  is,  the  density  is  constant  with  respect  to 
space.  Assuming  that  the  gas  velocity  at  the  gun  tube  entrance  is  zero,  it 
follows  that  the  gas  velocity  is  a  linear  function  of  distance.  Since  there 
is  a  mass  flow  into  the  gun  tube,  this  is  not  strictly  correct.  However, 
the  derivation  with  a  velocity  at  the  entrance  is  much  more  complicated,  and 
the  standard  Lagrange  distribution  is  used  as  a  first  approximation.  This 
leads  to 


m 

OP 

m 


P(X)  =  Pj^  +  (Pj^  -  Pj^g) 


Where  is  the  distance  from  the  tube  entrance  to  the  base  of  the 
projectile.  Integrating  over  the  tube  length. 


PrS  3M  . 


I  ^  3M 


where  p^  is  the  average  pressure  in  the  tube.  This  can  be  rearranged  as 


-  "rs  V‘“£i> 

^R  1  +  M,7(3M  .) 

4  pi 


Pr  >  Prs 


where  pj^  is  the  pressure  at  the  right  end  of  the  gun  tube  (base  of  the 
projectile),  is  the  mass  of  the  gas  in  the  tube,  Mpj  is  the  mass  of  the 

projectile,  and  pj^g  is  the  resistance  pressure.  The  latter  is  an  input 
parameter,  which  takes  into  account  the  shot  start  engraving  force  and 
frictional  forces  between  the  projectile  and  the  bore.  The  pressure  at  the 
gun  tube  entrance  is 


Pl  PrI^  2m  .1  Prs  2m  . 

\  pj  /  pj 


Pr  >  Pf 


The  above  two  equations  assume  that  the  projectile  has  started  to  move, 
creating  a  pressure  differential.  If  the  projectile  has  not  moved,  then 
region  4  is  treated  as  just  an  extension  of  region  3,  that  is. 


Given  the  pressure  on  the  projectile,  the  acceleration  equation  is 
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dV  . 

>*“i  -df^  ■  ‘fR 


'■rS*  *4  8o'«pJ 


P  >  P 
R  RS 


(A3) 


The  projectile  travel  Spj  is  then  given  by 


’^R  ^  ^RS 


dS  . 

1“1  -dl'  -  • 


(AA) 


The  final  two  differential  equations  are  exactly  analogous  to  region  3. 
That  is. 


dT'-vI  dT^vf 

H  H 


(A5) 


where  is  the  density  of  the  gas  in  region  A,  and 


^Pa  "^^a 

‘“I  di--i;-dJ-  " 


V,  -  bM^ 


m,,  >  0 
34 


(46) 


io  dt 


“3A  <  « 


where  is  the  mass  flux  into  region  A,  given  by 


Sa  '  S  ^a’^^®o  ^3^P3  ■  ^A^ 


P3  >  Pa 


(A7) 


%  \  /  2g^  p^(P^  -  P^) 


P3  <  Pa  • 
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The  speed  of  sound  in  region  4  is  given  by 


___foJP4__ 


IV.  VENT  OPTIONS 

The  model  described  above  has  a  simple  vent  option.  The  vent  area 
between  the  liquid  chamber  and  the  combuslon  chamber  is  fixed,  and  equals 
the  total  area  of  the  holes  drilled  in  the  piston.  Three  vent  options  are 
included  in  the  code,  where  the  constant  vent  area  is  a  special  case  of 
option  1.  For  each  option,  the  input  and  the  vent  area  calculation  is  in  a 
separate  subroutine,  and  other  options  can  be  added  easily. 

A.  VENT  1 

A  table  of  vent  areas  is  read  in  as  a  function  of  piston  travel. 

Linear  interpolation  is  used  to  determine  the  vent  area  for  any  given  piston 
travel.  The  case  of  constant  vent  area  is  a  special  case  of  this  option. 

B.  VENT  2 

A  standard  concept  in  Interior  ballistics  is  that  of  the  optimum 
(constant  pressure)  gun.  That  is,  for  a  given  gun  design,  we  want  to  know 
the  maximum  possible  performance  (muzzle  velocity)  for  a  given  limit  on  the 

9 

breech  pressure.  There  can  also  be  an  acceleration  limit  on  the  projectile 
(base  pressure).  For  a  lumped  parameter  solid  propellant  gun,  this  is 
fairly  straightforward.  These  models  usually  assume  a  Lagrange  pressure 
distribution  or  a  similar  fixed  pressure  distribution.  There  is  a  known 
relation  between  the  breech  pressure  and  the  pressure  at  the  base  of  the 
projectile.  Usually,  the  shot  start  pressure  is  chosen  as  equal  to  the 
desired  base  pressure,  so  the  projectile  does  not  move  until  the  desired 
pressure  is  reached.  The  burning  rate  for  the  propellant  is  adjusted  so  as 


to  reach  the  desired  maximum  breech  (base)  pressure  and  remain  at  this  level 
until  propellant  burnout. 

For  a  liquid  propellant  gun,  the  relation  between  the  breech  pressure 
and  the  base  pressure  Is  not  known  analytically,  because  of  the  added 
complication  of  the  piston  motion.  Also,  the  burning  rate  cannot  be 
controlled  directly,  since  this  depends  on  the  mass  flux  Into  the  combustion 
chamber.  What  can  be  controlled  Is  the  vent  area. 

Our  scheme  assumes  that  there  is  a  limit  both  on  the  breech  pressure 
and  the  base  pressure  (projectile  acceleration).  The  input  parameters  are 
AVMIN  ■  initial  vent  area  *  minimum  vent  area,  AVMAX  “  maximum  vent  area, 
and  ACS  =  maximum  acceleration  (kllo-g).  The  desired  base  pressure  Is 

Pjj  -  Mpj  (.098146)  ACS/A^  -  (49) 

where  the  number  .098146  Is  a  conversion  constant.  Then  we  define 

RATIO  =  (Pjj  -  Pr)/Pd  (50) 

and  let 

A^  (new)  -  A^  (old)  (1  +  AK*RATIO)  (51) 

where  AK  Is  a  proportionality  constant  (usually  about  one).  However,  the 
new  vent  area  Is  overridden  by  AVMIN  and  AXMAX. 

In  practice,  the  vent  area  stays  at  AVMIN  until  the  desired  base 
pressure  has  been  achieved.  The  projectile  then  starts  moving,  and  the  base 
pressure  falls.  The  vent  area  then  Increases  over  several  time  steps  to 
maintain  the  base  pressure.  There  Is  a  delay  In  Increasing  the  vent  area 
and  a  further  delay  until  the  pressure  Increase  In  the  combustion  chamber 
affects  the  pressure  distribution  In  the  gun  tube.  So  the  base  pressure 


falls  and  then  Increases  again.  AVMAX  is  normally  chosen  by  trial  and  error 
so  that  the  breech  pressure  does  not  exceed  a  desired  limit.  Once  AVMAX  is 
reached,  the  breech  and  base  pressures  again  fall  off. 

There  appears  to  be  no  way  to  implement  exactly  a  constant  pressure 
algorithm  for  the  regenerative  gun  without  major  modifications  in  the  piston 
motion  algorithm.  The  above  algorithm  does  give  an  approximation  to  the 
maximum  performance  possible  for  given  restrictions  on  the  breech  pressure 
and  the  base  pressure. 


C.  VENT  3 

Recent  regenerative  guns  have  been  designed  with  an  annular  piston  (see 
fig.  2).  The  piston  has  a  circular  hole  in  the  center  around  a  central  rod 
or  bolt.  The  fuel  is  stored  between  the  bolt  and  the  piston.  The  bolt 
remains  fixed  as  the  piston  moves.  The  motion  of  the  piston  opens  a  wider 
annular  gap,  and  fuel  is  forced  between  the  piston  and  the  bolt. 

We  assume  that  the  bolt  is  made  of  frustrums  of  cones.  Let  be  the 
distance  from  the  combustion  chamber  end  of  the  bolt  to  the  junction  between 
piece  1-1  and  1.  Then  the  volume  of  the  bolt  is 


TT  2  2 

5  ^’'i  “  *i-l^^’'^l^  ^i*^i-l  ’^i-1  ^  * 


Again  we  assume  that  the  piston  has  zero  thickness.  Then  the  volume  of  the 
liquid  is 


Vj  .  ,r^  (x^  -  Xj)  - 


r  ir(x.  -  X  )  (r 

1.2  ^  ^ 


I  (^i  +  r^r^_j  +  )J 


y  V  i  "II  rv  i-_k  aIV.  r  *  •_ 


*  •  *  «  * 


where  rj^  is  the  radius  of  the  liquid  chamber  and  W-1  is  the  number  of  pieces 
in  the  bolt. 

The  required  information  to  determine  the  geometry  is  the  and  the  r^ 
values.  In  practice,  we  often  want  to  specify  the  total  liquid  volume  and 
the  vent  areas  for  parametric  studies.  So  the  actual  input  used  is  a  table 
of  distances  and  vent  areas,  as  well  as  the  volume  and  area  of  the  liquid 
chamber,  and  the  area  of  the  hole  in  the  piston.  The  radii  and  then  the 
volume  of  the  chamber  can  be  computed.  By  equation  (53),  the  volume  is 
directly  proportional  to  the  distances  Xj^.  So  the  x^  are  scaled  to  give  the 
desired  Vj.  If  the  dimensions  of  an  actual  gun  are  used,  where  the  input 
liquid  volume  matches  the  input  dimensions,  no  adjustment  is  made. 

There  are  several  changes  necessary  in  the  previous  governing 
equations.  Let  the  piston  travel  satisfy  the  relation  <  Sp^  <  x^. 

The  radius  of  the  rod  at  this  point  can  be  found  by  linear  interpolation 


■■  ■  ■'i-i  + 


(S  -  X.  ,) 

:_E£ _ ili.  (r  _  r  ) 

(Xi  -  x^_p  ^’'i  Vl'*  • 


Taking  the  derivative. 


dt  ps  (x^  -  x^_^) 


Now  the  volume  of  the  liquid  reservoir  is  given  by 


''l  “  '^lo  “  ^Ss  ■  I  ^^s  ■  ^  ’^i-l  ^ 


+  1  5  (X 

j-2  J 


"j  "j-1 


or  taking  derivatives 


odSdS- 
.2  PS  ^  IT  ps  ,  2  .  ,  2v 

L  "dt~  3  “dt  ^  ’^i-l  ’^i-l  ^ 


m 


"  5  (^s  -  *i-l^  <2r  .  r^_^  ^I)  . 


Rearranging  this. 


'll  ^  -  -  *1  *  5  pps  '  'l-l  * 

■"  <%s  ■  -l-l’  +  tj.j)  . 


Similarly, 


dr  “  '^ps  ^3 


PS  ^-1  * 


*  '%s  '  ’‘l-l>  ' 


) 

i-r  dtj 


The  piston  acceleration  is  given  by 


dV  g 

^2]  “dt“  *  m“  ^*^3  ^^3  '  ~  **1  ^^1  ~ 

ps 


where  is  the  area  of  the  hole  in  the  piston.  The  acceleration  on  the 
piston  no  longer  depends  on  the  vent  area.  The  vent  area  Itself  is  given  by 
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Ay^Spg)  *  Ajj  -  Ajj(Spg) 


(61) 


This  is  used  to  compute  the  mass  flux  other  governing  equations 

are  unchanged. 


V.  DROPLET  BURNING  OPTIONS 


So  far  we  have  assumed  that  the  liquid  combusts  instantaneously  upon 
entering  the  combustion  chamber.  There  is  some  evidence,  however,  that 
liquid  accumulates  in  the  combustion  chamber.  To  see  what  effect  this  may 
have  on  the  gun  performance,  simple  approximations  for  the  formation  and 
combustion  of  droplets  in  the  combustion  chamber  are  considered. 

The  behavior  of  the  liquid  as  it  enters  the  combustion  chamber  is  very 
complicated.  We  will  assume  that  the  liquid  instantaneously  forms  droplets 
of  a  fixed  size.  The  droplets  will  then  combust  at  some  given  rate.  The 
initial  size  of  the  droplets  is  an  input  parameter. 

In  the  above  equations,  we  included  the  possibility  that  gas  could  flow 
from  the  gun  tube  back  into  the  combustion  chamber.  In  practice,  this 
occurs  only  rarely,  and  the  amounts  Involved  are  very  small.  Since  the 
equations  below  are  already  very  complicated,  an  additional  assumption  is 
made  that  the  mass  flux  from  the  combustion  chamber  to  the  gun  tube  is 
positive.  This  eliminates  some  of  the  algebra,  and  has  no  noticeable  effect 
on  the  results. 


A.  DROP  1 

The  liquid  combusts  Instantaneously  as  it  enters  the  combustion 
chamber,  releasing  all  of  its  chemical  energy.  This  is  the  option  derived 
above . 
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%  _%  . 


B.  DROP  2 


Gough^  has  implemented  a  simple  droplet  combustion  scheme.  The 

droplets  have  a  fixed  size  at  all  times,  and  combustion  simply  reduces  the 

number  of  drops.  The  effect  of  droplet  accumulation  may  be  considered, 
without  implementing  a  whole  distribution  of  droplet  sizes.  We  will  derive 
the  appropriate  equations  for  this  scheme  in  our  lumped  parameter  model. 

Let  be  the  liquid  mass  in  region  3,  liquid  volume, 

be  the  gas  mass,  and  Vg3  be  the  gas  volume.  Then 

'“u  - 


and 


^G3  “  “g3'^'^G3 


(63) 


where  is  the  density  of  the  liquid  and  is  the  density  of  the  gas. 
Define  the  porosity 

'3  ■  ''C3^''3  • 


Then 


^3 


M3/V3 


(65) 


-  (1  -  £3)  Pj^3  +  C3  pQ3  . 


Following  eq.  (17),  the  speed  of  sound  in  the  liquid  is 


c 


L3 


(66) 
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and  following  eq.  (27)  the  speed  of  sound  in  Che  gas  Is 


■G3 


/nzi!3 

"03  “  '  ^  '>G3> 


(67) 


The  liquid  Is  still  considered  Co  be  Isothermal,  so  the  Internal  energy  does 
not. -change.  The  heat  transfer  to  the  liquid  is  Ignored  because  of  the 
complexity  of  Including  this  effect.  The  enthalpy  of  the  liquid  Is 

where  e^  Is  the  chemical  energy  of  the  liquid.  Let  m2  be  the  rate  at  which 
the  liquid  in  region  3  Is  combusting  and  forming  the  final  gas  products. 

The  thirteen  equations  derived  In  Section  III  are  the  same  except  for 
the  energy  equations  [8]  and  [13].  To  derive  the  new  equations,  we  need 
some  preliminary  results.  Consider  the  mass  conservation  equations 


dt 


d 

dt 


[Pj^3(l  -  €2)  V2] 


\3  . 

m”  “34  ■  “3 


(69) 


and 


dM 


G3 

dt” 


dt  ^3  ^3^ 


G3  . 
m.  -  - —  m. . 


(70) 


Expanding  the  equations, 


"dt  "  ^3^  ^^3  ~  ^L3'^3  dt  ^3 


\3  . 


“l3  "  M~  “34  "  “3 


‘*"g3  „  ,  V  ‘^!3.  .  f!3 

"dt~  ®3'^3  ^G3  ^^3  "dt  ^G3  3  dt 


.  G3  . 

*  ”3  ■  "M3  “34  * 


Combining  the  two  equations, 


Now  consider  the  pressure  equations.  Assume  that  the  gas  and  liquid 
pressures  In  the  region  are  Identical.  The  liquid  pressure  equation 
following  eq.  (18),  is 


‘**’3  ‘^L3  ‘*^L3 


dt  g^  dt 


and  the  gas  phase  equation,  following  eq.  (33),  is 


m 


*^^3  *^03  *^^03  .  “3^*^L3  "  "g3^^^  ”  . 


dt  g^  dt 


''g3  -  “G3 


Note  the  gas  phase  equation  Is  no  longer  affected  directly  by  the  liquid 
Injection  Into  the  chamber,  but  by  the  rate  at  which  the  liquid  Is  turned 
Into,  gas.  Now  we  define  the  speed  of  sound  In  the  mixture  by  the  equation 


1  ^3  ~  ^3^ 

2  “  2  2" 


^3*^3  ‘’g3‘^G3  ^L3‘^L3 


i_  r _ 1 

^3  e^/(pQj  0^3”)  + 


+  (1  -  e3)/(PL3  ^^3 


Combining  eq.  (74)  and  (75),  and  using  eq.  (76),  the  result  Is 


‘"’3.V3Lf!3,S3 

dt  V3  dt 


S.  1  .  !c3  ^  fo 
^G3  ^3 


o  ^^3  "■  **G3^^^  ’  *34 


Cg3  (1  -  bp^^) 


Now  consider  region  4.  The  mass  conservation  equations  are 


**^4  \3  .  . 

"dt"  “  “34  4 


G3  . 

”  "34  ^  "4 


or  expanding  as  before 


'dt"  ■  '4*  ''4  '  "iA  '^4  dt' 


*  „  ,  \3  . 

*  ‘>L4  <*  -  '4>  -dJ  ■  S;-  "34  -  "4 


‘*^04  ‘^‘=4  ‘*'^4  ”g3 

-di-  Pg4  ^  -dl  ^  ^G4  S  -dl  -  "34  "4  • 


Again  we  assume  that  the  gas  and  liquid  pressures  In  the  region  are 
Identical.  Then  the  pressure  equations  are 


dt  ”  g  dt 
o 


'g4  '‘‘>04  .  ‘4<\4  -  '’g4>'’'  - 


G3  34'  G3 


Note  that  the  mixture  In  region  3  Is  assumed  to  be  homogeneous,  and  the 
fluid  that  flows  Into  region  4  has  the  same  porosity  as  the  fluid  In  region 
3.  In  general,  region  4  has  a  different  porosity.  So  the  gas  phase 


pressure  equation  for  region  4  has  to  Include  the  effect  of  that  portion  of 
the  mixture  in  region  3  that  Is  gas.  Combining  the  equations  before 
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dP4  P4  c 

“31  dT-r-v- 

o  4 


^^4  ®4 
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(85) 


.  _“34. 

^04^3 


^G4  ,  ®o  ”g3  ^^L3  “  ^4^^^  " 


M„,  +  Mj^3  - 


G3 


■'L4 


‘^GA  -  '’Pg4> 


To  complete  the  system,  equations  are  required  to  determine  the  liquid 
density  and  mass  In  the  regions.  The  corresponding  gas  quantities  can  be 
easily  derived.  The  liquid  pressure  equations  (74)  and  (83)  are  rewritten 
as 


[14] 


‘*^13  So  ‘‘P3 

”dt“  ''"I  'dt 

L3 


and 


[16] 


(86) 


(87) 


The  mass  conservation  equations  for  the  liquid  are  written  as 


[15] 


“dt 


\3  . 


®13  ■  "S'  “34  ”  “3 


and 


(88) 
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(17] 


(89) 


!!l3  . 

”dt  “34  ”  “4  • 


To  close  the  system,  we  need  information  on  the  rate  at  which  the 
liquid  droplets  are  combusting.  The  rate  of  surface  regression  is  assumed 
to  be  of  the  form 

rate  of  surface  regression  =  A  p  .  (90) 

Liquid  propellant  burning  rates  have  been  measured  by  McBratney.^^  The 
rate  of  burning 


*3  - 


■^L3 


S  A  p 


(91) 


where  S  is  the  total  surface  area  of  the  droplets  in  a  region.  All  the 
droplets  are  assumed  to  have  a  constant  diameter  d. 

The  surface  area  of  a  single  drop  is  given  by 


S 


D 


ird 


(92) 


and  the  volume  of  a  single  drop  by 


V 


D 


ir 

6 


(93) 


^^W.F.  McBratney,  "Windowed  Chamber  Investigation  of  the  Burning  Rate  of 
Liquid  Monopropellants  for  Guns,"  ARBRL-MR-03018 ,  April  1980. 

^^W.F.  McBratney,  "Burning  Rate  Data,  LPG  1845,"  ARBRL-MR-03128 ,  August 
1981. 


38 


If  Np  is  the  total  number  of  drops  in  region  3,  then 


"d  ■  ■ 


So  the  total  surface  area  in  the  region 


S  -  6V^3/d 


and  eq.  (91)  can  be  written  as 


“3  "  ^L3  '^L3  d  ^^3 


*3  ^3  d  ^^3 


For  region  4,  the  pressure  is  no  longer  constant  over  the  region,  but 
follows  the  Lagrange  pressure  distribution.  But  this  would  be  quite 
complicated  to  keep  track  of,  and  very  little  liquid  goes  into  region  4  for 
most  problems.  So  the  combustion  rate  is  assumed  to  depend  on  the  average 
pressure  P^,  and 


*4  •  "la  I  ‘Pi' 


-  -  A*,  -  ■'.V  ••  >  ■ 


c 


DROP  3 


The  above  representation  assumes  that  all  the  droplets  remain  the  same 
size  at  all  locations  and  times.  In  practice,  some  distribution  of  droplets 
will  exist  as  a  function  of  space  and  time.  In  this  section,  we  present  a 
slight  generalization  of  the  above  model. 

.  The  drops  still  form  Instantaneously  as  the  fluid  enters  the  combustion 
chamber,  and  all  the  droplets  are  initially  the  same  size.  However,  the 
droplets  will  shrink  as  they  burn.  They  will  be  represented  by  an 
approximation  to  a  probability  distribution.  One  advantage  is  that 
different  probability  distributions  can  be  easily  implemented  using  the 
formalism  developed  below.  So  alternate  theories  of  droplet  formation  and 
combustion  can  easily  be  added. 

A  standard  procedure  is  to  define  a  probability  distribution  function 
f(r)  *  the  probable  number  of  droplets  with  radius  between  r  and  r  +  dr.^^ 
For  our  purposes  it  is  more  convenient  to  use  the  mass  of  the  droplets 
rather  than  the  number  of  droplets.  So  define  M(r)  =  the  mass  of  the 
droplets  in  region  3  with  diameter  between  r  and  r  +  dr.  Since  this  is  a 
lumped  parameter  region,  there  is  no  space  dependency.  Then  let  m(r)  be  the 
rate  of  consumption  of  drops  with  radius  r,  or 


m(r)  =  M(r)  ^  Ap^®  . 


(99) 


Then  we  have  the  partial  differential  equation 


3M(r) 

3t 


A 

-  -m7  “34 


m(r) 


W(r2  3r 
3r  3t 


(100) 


^^F.A.  Williams,  "Progress  in  Spray-Combustion  Analysis,"  in  The  Eighth 
International  Combustion  Symposium,  Williams  and  Wilkins  Co.,  Baltimore,  MD, 


The  first  term  is  the  loss  term  to  region  4,  the  second  term  is  the 
loss  term  to  the  gas  phase,  and  the  last  term  is  the  loss  to  smaller 
drops.  There  is  also  an  input  boundary  condition  for  r  =  initial  radius  of 
the  droplets  and  an  output  boundary  condition  at  r  =  0. 

Rather  then  solve  the  above  partial  differential  equation,  the 
distribution  is  discretized  to  obtain  a  set  of  ordinary  differential 
equations.  To  be  consistent  with  the  previous  notation,  the  diameter  rather 
than  the  radius  of  the  droplets  will  be  used.  More  or  less  arbitrarily,  the 
distribution  is  split  into  ten  classes.  Let  *  the  initial  diameter  of 

the  droplets  and  d^^^^  =  be  the  difference  in  diameter  of  the  ten 

droplet  classes.  Then  define  “  the  mass  of  the  droplets  with 

diameter  between  -  (i)dinc  ‘^init  ”  ^^"'■^^‘^inc*  ‘^i  ° 

(1.05  -  0*li)diriit  the  working  diameter  for  each  class  of  droplets.  That 
is,  for  purposes  of  obtaining  burning  rates,  assume  that  this  is  the 
diameter  of  all  droplets  in  the  class.  Then 


msCi)  -  ^^3 


(101) 


is  the  mass  burning  rate  for  each  class,  and 


10 

m  =  I  m„(i) 
i=l 


(102) 


is  the  total  mass  burning  rate. 

To  derive  the  governing  ordinary  differential  equations,  assume  that  in 
any  class,  the  mass  is  evenly  distributed  among  the  possible  diameters. 

This  will  contradict  the  assumption  above  that  all  the  droplets  have  an 
average  diameter.  But  in  the  limit  of  infinitely  many  classes,  it  will  be 
equivalent.  Looking  at  the  first  class. 
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(103) 


M, 


or  the  change  in  mass  of  class  1  -  Influx  Into  region  3  -  efflux  to  region  4 
-  loss  to  the  gas  phase  -  loss  to  the  next  smaller  class.  Since  we  assume 
that-  the  mass  is  equally  distributed  over  the  possible  diameters, 

Ad/dinc  of  the  mass  will  shrink  into  the  next  smaller  class.  Taking  the 
limit 


I18J 


dMi3(l)  _  M,3(U  .  .  , 

"dt  ■  “13  '  "kT'  34  ■  57  ^*>’3  • 

i  inc 


(104) 


Note  that  the  rate  of  change  of  the  diameter  ■  twice  the  surface 
regression  rate  -  2Ap2  .  Since  we  assume  that  all  the  droplets  initially 
start  in  class  one,  all  the  Influx  from  region  1  goes  into  class  1.  For  the 
other  classes,  the  equations  are 


[19-27] 


dM,_3(i) 

dl 


\3<“ 

H„ 
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(105) 


+ 


M^3(i  -  1)  -  Mj^3(1) 

d, 

Inc 


2  ...10. 


All  the  earlier  equations  [1]  to  [17]  from  the  above  section  are  the 
same,  except  for  the  new  definition  of  m^  and  the  corresponding  definition 
of  m^.  One  of  the  equations  in  (105)  is  in  fact  redundant,  since  * 
Mj^sd)  +  ...  M^3(10).  The  extra  equation  is  included  for  the  sake  of 
simplicity. 

The  equations  for  region  4  are  similar,  except  that  there  is  no  source 
term  from  region  1.  The  loss  terms  for  region  3  are  the  source  terms  for 


r 


region  4.  That  is, 


[28-37]  — g-  -  m^d) 
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VI.  MASS  AND  ENERGY  BALANCE 


As  a  check  on  the  equations  derived,  both  mass  and  energy  should  be 
conserved.  The  mass  balance  Is  straightforward. 


M.p  -  Mj  +  M3  + 


where  Mj  Is  the  total  mass  of  the  propellant  and  primer.  This  should  be 
constant  throughout  the  integration. 


The  energy  balance  is  more  complicated.  The  energy  loss  through 
friction  is  quite  small,  and  is  ignored.  The  code  presently  does  not 
include  heat  loss  or  air  resistance.  The  liquid  is  considered  to  be  an 
isothermal  fluid,  and  its  internal  energy  e  is  Just  the  chemical  energy  of 
the  propellant.  The  total  energy  of  the  liquid  is 


El  *  ®i^Mi  +  Mj^^  +  Mj^)  . 


The  internal  energy  of  the  gas  in  region  3  is  given  by 


'3  ■  'v  ■'3 


and  the  total  gas  energy  is 


m 

IfT" 


h  ■  'v  ■'3  “g3  • 


The  average  Internal  energy  of  the  gas  In  region  4  is  given  by 


e  -  c  T,  dx 

4  X  J  o  V  4 
R 


Clio; 


(111) 


e  -  i  /  ^  - 7 - -pr—  dx 


(112) 


by  the  Noble-Abel  equation.  Using  eq.  (38)  for  the  pressure  function,  this 
can  be  Integrated  to  yield 


(1  -  bp^) 

®4  “  P4 


(113) 


^4  ”  ®4”g4  . 


The  kinetic  energy  of  the  piston  is  given  by 


ps  ps  ps  o 


and  the  kinetic  energy  of  the  projectile  by 


EK  ^  -  0.5M_,  V_//g_. 


(114) 


(115) 


(116) 


By  our  assumptions,  only  the  fluid  In  region  4  Is  moving.  The  kinetic 
energy  of  the  gas  Is 


^^G4  “  /o^  0-^  Pg4  \  * 


We  Ignore  drag  and  assume  that  the  liquid  has  the  same  velocity  as  the 
gas.  Then  the  kinetic  energy  of  the  liquid  Is 


The  Lagrange  pressure  distribution  Implies  a  linear  velocity  profile  with 
the  gas  velocity  at  the  tube  entrance  equal  to  zero.  This  Implies  that 


*  -  Vpj  x/x^  . 


Substituting  and  Integrating 


®''g4  ■  "-S  "g4 


EKl4  ■  0.5  Vpj'  x^/(3g„)  . 


Then  the  total  energy  of  the  system  Is 


E_  -  E,  +  E-  +  E,  +  EK„„  +  EK  ,  +  EK_,  +  EK, , . 
T  1  3  4  ps  pj  G4  L4 
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This  should  be  constant  throughout  the  integration 


VII.  NUMERICAL  METHOD 

The  ordinary  differential  equations  derived  above  are  solved  using  the 
code  EPISODE.^®  This  is  a  robust  and  efficient  code  for  the  solution  of 
ordinary  differential  equations.  EPISODE  proper  consists  of  the  subroutines 
DRIVE  to  SING  (see  Appendix  A). 

Only  two  changes  have  been  made  in  EPISODE  proper.  The  first  concerns 
the  error  control.  Originally  this  could  be  either  a  relative  or  an 
absolute  error  control.  A  relative  error  control  can  be  wasteful  if  some  of 
the  quantities  being  integrated  become  negligibly  small.  However,  an 
absolute  error  control  is  inaccurate  if  the  quantities  integrated  vary 
widely  in  magnitude. 

A  solution  to  this  problem  was  developed  during  work  on  integrating 

chemical  kinetics  networks,  where  the  concentrations  of  species  differ  by 

1  q 

many  orders  of  magnitude.  This  is  a  semi-relative  error  control.  If  the 
quantities  are  above  some  cutoff  value  (SREC)  a  relative  error  control  is 
used.  If  they  are  below  SREC,  an  absolute  error  control  is  used.  While  an 
absolute  error  control  would  be  adequate  for  the  simpler  equations  now  being 
considered,  the  semi-relative  error  control  has  been  left  in  the  code.  SREC 
has  normally  been  chosen  as  10~®. 

The  other  change  is  at  the  end  of  the  routine  TSTEP,  which  actually 
takes  the  integration  steps.  Several  options  have  been  included  here. 

First,  diagnostic  printouts  of  the  quantities  of  interest  may  be  made  after 


^®A.C.  Hindmarsh  and  G.D.  Byrne,  "EPISODE;  An  Effective  Package  for  the 
Integration  of  Systems  of  Ordinary  Differential  Equations,"  UCID-30112- 
Rev.l,  Lawrence  Livermore  Laboratory,  1977. 

^^T.P.  Coffee,  "A  Computer  Code  for  the  Solution  of  the  Equations  Governing 
a  Laminar,  Premixed,  One-Dimensional  Flame,"  ARBRL-MR-03165 ,  April  1982. 
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each  time  step  (KWRITE*!).  Second,  if  the  optimum  gun  option  is  being  used 
(VENT  2),  the  vent  area  is  updated  at  the  end  of  each  time  step.  Finally, 
some  information  about  regions  1  and  4  must  be  included.  Eq.  [4]  for  the 
liquid  density  becomes  singular  as  Vj  approaches  zero.  So  when  becomes 
some  fraction  VPER  of  the  initial  volume,  the  region  is  closed.  That  is, 
region  1  is  ignored  in  the  integration  from  then  on.  This  will  lead  to  a 
small  error,  due  to  the  liquid  propellant  left  in  region  1  when  the  region 
is  closed.  However,  VPER  is  an  input  parameter,  and  by  trial  and  error  can 
be  chosen  small  enough  that  the  error  is  negligible.  If  ever  becomes 
negative,  the  step  is  rejected,  and  the  integration  starts  again  with  a 
smaller  time  step.  There  is  a  similar  problem  with  region  4,  since  the 
projectile  may  start  at  the  throat,  and  will  be  zero.  So  region  4  starts 

out  as  closed.  When  the  projectile  is  at  least  .01  cm.  down  the  gun  tube, 
the  region  is  opened.  Before  this,  region  4  is  considered  as  an  extension 
of  region  3. 

There  are  actually  two  integration  routines  included  in  EPISODE.  The 
first  is  the  Adams  method,  suitable  for  non-stiff  problems  (METH=1).  The 
second  is  the  backward  differentiation  method,  suitable  for  stiff  problems 
(METH=2).  We  have  normally  used  METH*2.  However,  the  governing  equations 
are  not  very  stiff,  and  METH=1  can  also  be  used.  The  run  times  are  about 
equivalent  for  the  two  methods. 

There  are  also  several  methods  for  the  solution  of  the  nonlinear  system 
of  equations  required  at  each  step  of  the  integration.  We  have  used  the 
Newton  method  with  an  internally  computed  finite  difference  approximation  to 
the  Jacobian.  So  the  only  subroutine  that  needs  to  be  supplied  by  the  user 
is  DIFFUN.  Given  a  set  of  values  for  the  unknowns,  DIFFUN  computes  the  time 
derivatives  of  the  unknowns. 

Subroutine  DIFFUN  in  this  code  only  calls  one  of  three  subroutines 
(FDROPl,  FDR0P2,  or  FDR0P3).  These  routines  correspond  to  the  instantaneous 
burning  option,  the  fixed  drop  size  option,  and  the  shrinking  drop  size 
option.  Since  the  governing  equations  are  different,  it  is  easier  to  set 
these  up  as  separate  subroutines. 


The  rest  of  the  routines  have  to  do  primarily  with  input  and  output. 

The  main  routine  RLGD  read  in  the  input  data,  sets  the  initial  conditions, 
calls  the  integrator,  and  writes  the  main  output  file. 

To  make  it  easier  to  add  or  change  options,  some  of  the  input  is  in 
separate  subroutines.  For  instance,  subroutines  VENTl ,  VENT2,  and  VENT3 
read  in  the  required  data  for  the  three  vent  options.  These  routines  also 
compute  the  vent  area  for  a  given  piston  travel.  So  if  a  new  vent  option  is 
desired,  it  is  only  necessary  to  add  a  new  subroutine  and  change  the  few 
lines  that  call  the  VENT  subroutines. 

Similarly,  DR0P2  and  DROPS  read  in  the  input  parameters  for  the  droplet 
options.  PISRES  reads  in  the  piston  resistance  as  a  function  of  piston 
travel.  DISl  and  DIS2  read  in  the  discharge  coefficients.  These  can  be 
either  a  function  of  piston  travel  or  a  function  of  liquid  pressure. 

PROJRES  reads  in  the  projectile  resistance  pressure  as  a  function  of 
projectile  travel.  Subroutine  PRIM2  is  a  simple  primer  option.  That  is, 
the  primer  may  be  included  as  droplets  in  region  3  to  follow  the  initial 
behavior  of  the  system  in  more  detail. 

Subroutines  CAPTION  and  OUT  control  the  output  files.  Subroutine 
OUTGRA  creates  a  graphics  file,  so  that  any  of  the  quantities  computed  may 
later  be  graphed. 


VIII.  NUMERICAL  COMPARISONS 

As  a  check,  we  will  compare  our  code  with  the  numerical  code  developed 
by  Gough. ^  The  sample  problem  used  is  discussed  in  ref.  [6].  This  is  a  25 
mm.  gun  with  a  constant  vent  area.  The  input  data  required  is  reproduced  in 
Table  1. 


48 


Table  1.  Baseline  Input  Data  for  Test  Problems. 


fuel  side  area  of  piston 

13, 

.00 

2 

cm 

combustion  side  area  of  piston 

15. 

.48 

2 

cm^ 

gun  tube  area 

4, 

.91 

2 

cm 

initial  fuel  reservoir  volume 

96, 

.00 

3 

cm"' 

initial  combustion  chamber  volume 

35, 

.40 

3 

cm-* 

initial  projectile  offset  in  tube 

1, 

.27 

cm 

projectile  travel 

213. 

.36 

cm 

projectile  weight 

194, 

.40 

gm 

piston  weight 

480, 

.82 

gm 

bore  discharge  coefficient 

.95 

piston  discharge  coefficient 

.58 

shot  start  pressure 

36, 

.00 

MPa 

liquid  density  (1  atm.) 

1 

.23 

gm/cm 

bulk  modulus 

2413 

.10 

MPa 

liquid  chemical  energy 

3330 

.90 

joules/gm 

ratio  of  specific  heats 

1 

.26 

molecular  weight  of  combustion  gases 

19 

.01 

gm/mole 

covolume 

1 

.26 

cm^/gm 

initial  liquid  pressure 

.10 

MPa 

initial  gas  pressure 

81 

.00 

MPa 

The  answer  given  here  for  the  Gough  code  are  slightly  different  than 

those  reported  In  ref.  6.  The  Gough  code  controls  the  Integration  time  step 

size  based  on  the  grid  In  the  gun  tube.  The  time  step  must  be  small  enough 

20 

not  to  violate  the  Courant-Frledrich-Levy  condition.  The  maximum  number 
of  grid  points  and  the  minimum  grid  spacing  in  the  gun  tube  are  input 
conditions.  A  smaller  grid  size  than  was  used  in  ref.  6  was  required  to 
obtain  convergence. 


^®P.J.  Roache,  Computational  Fluid  Dynamics.  Hermosa  Publishers,  1972. 
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The  results  for  these  computations  are  given  in  Table  2.  The  agreement 
between  the  two  codes  is  not  very  good. 


Table  2.  Comparison  -  Baseline  Input  Parameters. 


- 

Gough 

Coffee 

Muzzle  velocity  (m/s) 

1097 

1161 

Max.  liquid  pressure  (MPa) 

289 

351 

Max.  combustion  pressure  (MPa) 

231 

277 

Time  until  projectile  exit  (ms) 

3.7 

3.2 

Computer  time  (s) 

20.6 

1.2 

This  can  be  partially  explained  by  considering  the  way  the  Gough  code 
computes  the  mass  flux  between  the  liquid  chamber  and  the  combustion 
chamber.  If  the  combustion  chamber  pressure  is  higher,  gas  can  flow  into 
the  liquid  chamber.  Since  the  liquid  chamber  is  one  phase,  the  code 
converts  this  into  liquid,  which  defies  conventional  physics.  The  new  code 
normally  does  not  allow  flow  into  the  liquid  chamber.  However,  it  can  be 
run  with  this  option.  The  results  are  in  Table  3. 

Table  3.  Comparison  -  Back  Flow  into  Liquid  Chamber. 


Gough 

Coffee 

Muzzle  velocity  (m/s) 

1097 

1140 

Max.  liquid  pressure  (MPa) 

289 

316 

Max.  combustion  pressure  (MPa) 

231 

251 

Time  until  projectile  exit  (ms) 

3.7 

3.4 

Computer  time  (s) 

20.6 

0.7 

The  agreement  Is  better,  but  there  are  still  noticeable  differences. 
Another  possible  cause  for  error  Is  the  time  step  control  In  the  Gough  code, 
which  is  based  solely  on  the  conditions  in  the  gun  tube.  However,  in  the 
present  problem  we  are  starting  with  a  large  presrure  difference  between  the 
combustion  chamber  and  the  liquid  reservoir.  The  time  step  must  be  reduced 
to  resolve  the  interaction  between  these  two  chambers  and  the  piston 
movement . 

To  check,  this  hypothesis,  the  initial  pressure  conditions  were  changed 
to  initial  liquid  pressure  “  9  MPa  and  initial  gas  pressure  ■*  10  MPa.  This 
is  less  of  a  strain  on  the  error  control.  As  Table  3  shows,  the  two  codes 
now  give  almost  identical  answers. 


Table  4.  Comparison  -  Initial  Gas  Pressure  =  10  MPa. 


Gough 

Coffee 

Muzzle  velocity  (m/s) 

1114 

1133 

Max.  liquid  pressure  (MPa) 

311 

310 

Max.  combustion  pressure  (MPa) 

247 

247 

Time  until  projectile  exit  (ms) 

4.2 

4.2 

Computer  time  (s) 

23.5 

2.5 

The  only  noticeable  difference  is  in  the  muzzle  velocity.  The 
difference  is  in  the  expected  direction.  That  is,  the  new  code  assumes  a 
Lagrange  pressure  distribution  in  the  gun  tube.  So  any  changes  in  pressure 
in  the  combustion  chamber  Instantaneously  effect  the  projectile,  no  matter 
how  far  down  the  gun  tube  it  is  located.  The  Gough  code  uses  partial 
differential  equations  to  describe  the  gun  tube,  and  actually  follows  the 
pressure  waves  as  they  propagate.  So  there  is  a  delay  between  the  pressure 
rise  in  the  combustion  chamber  and  the  corresponding  pressure  rise  at  the 
base  of  the  projectile.  This  results  in  a  slightly  lower  muzzle  velocity. 


Fig.  3  compares  the  combustion  chamber  pressures  as  a  function  of  time, 
and  fig.  4  compares  the  projectile  velocities.  These  curves  confirm  the 
above  analysis.  The  agreement  is  excellent  at  early  times.  But  as  the 
projectile  moves  down  the  gun  tube,  differences  appear  and  become  greater 
with  time. 

It  is  useful  to  know  if  these  differences  can  become  important.  So  we 
took. the  previous  problem  and  increased  the  liquid  reservoir  volume  to  300 
cm  and  the  projectile  travel  to  1000  cm.  Results  are  given  in  Table  5  and 
in  figs.  5  and  6.  There  is  in  fact  a  fairly  large  difference  in  muzzle 
velocities . 


Table  5.  Comparison  -  Long  Gun  Tube. 


Gough 

Coffee 

Muzzle  velocity  (m/s) 

1556 

1685 

Max.  liquid  pressure  (MPa) 

264 

261 

Max.  combustion  pressure  (MPa) 

205 

203 

Time  until  projectile  exit  (ms) 

10.0 

10.2 

Computer  time  (s) 

32.8 

2.4 

In  conclusion,  the  Gough  code  may  be  inaccurate  if  the  mechanisms 
controlling  the  time  step  occur  in  the  lumped  parameter  regions.  The  new 
code,  however,  will  overpredict  the  muzzle  velocity  for  a  long  gun  tube.  It 
may  also  be  inaccurate  for  a  high  performance  (large  muzzle  velocity)  gun, 
if  the  projectile  moves  fast  enough  that  the  pressure  increase  in  the  throat 
of  the  gun  tube  cannot  catch  the  projectile.  The  new  code  runs  about  10  to 
15  times  faster  than  the  Gough  code.  Several  other  comparisons  have  been 
made,  and  the  above  conclusions  still  hold. 
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Figure  3.  Baseline  25  ram  Gun  -  Instantaneous  Combustion  -  Initial 
Chamber  Pressure  =  10  MPa.  -  Chamber  Pressures. 

Coffee  model  (line);  Gough  model  (dot) 


Baseline  25  mm  Gun  -  Instantaneous  Combustion  -  Initial 
Chamber  Pressure  =  10  MPa.  -  Projectile  Velocities. 
Coffee  model  (line);  Gough  model  (dot) 


Long  25  mm  Gun  -  Instantaneous  Combustion  -  Initial 
Chamber  Pressure  =  10  MPa.  -  Chamber  Pressures. 
Coffee  model  (line);  Gough  model  (dot) 


Long  25  mm  Gun  -  Instantaneous  Combustion  -  Initial 
Chamber  Pressure  =  10  MPa.  -  Projectile  Velocities. 
Coffee  model  (line);  Gough  model  (dot) 
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APPENDIX  A 
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A  listing  of  the  computer  code  follows.  Routines  RLGD  to  PEDERV  have 
been  written  for  the  regenerative  liquid  gun  application.  Routines  DRIVE  to 
SING  comprise  the  code  EDISODE. 
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APPENDIX  B 


Below  Is  a  listing  of  the  job  stream  for  the  above  test  problem 
(initial  gas  pressure  «  10  MPa),  Following  is  a  listing  of  the  output. 
Since  there  are  many  variables  to  keep  track  of,  the  output  is  on  multiple 
files.  Only  the  ones  of  Interest  for  a  particular  problem  need  to  be 
printed  out.  Additional  files  contain  the  values  of  interest  for  a  problem 
with  a  droplet  option. 
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GLOSSARY 


Cross  section  area  of  the  liquid  reservoir,  cm  . 

2 

Cross  section  area  of  the  combustion  chamber,  cm  . 

2 

Cross  section  area  of  the  gun  tube,  cm  . 

2 

Vent  area,  cm  . 

O 

Covolume,  cm-^/gm. 

Speed  of  sound  In  the  liquid  in  the  reservoir,  cm/s. 

Speed  of  sound  In  the  liquid  in  the  combustion  chamber,  cm/ 

Speed  of  sound  In  the  liquid  In  the  gun  tube,  cm/s. 

Speed  of  sound  in  the  gas  In  the  combustion  chamber,  cm/s. 

Speed  of  sound  in  the  gas  In  the  gun  tube,  cm/s. 

Average  speed  of  sound  in  the  combustion  chamber,  cm/s. 

Average  speed  of  sound  In  the  gun  tube,  cm/s. 

Specific  heat  at  constant  pressure,  joules/gm-K. 

Specific  heat  at  constant  volume,  Joules/gm>K. 

Discharge  coefficient  for  the  piston. 

Discharge  coefficient  for  the  gun  tube. 

Chemical  energy  of  the  liquid,  Joules/gm. 


Internal  energy  of  the  gas  In  the  combustion  chamber,  joules/gm. 
Internal  energy  of  the  gas  in  the  gun  tube,  joules/gm. 

Total  internal  energy  of  the  liquid.  Joules. 

Total  Internal  energy  in  the  gas  in  the  chamber,  joules. 

Total  Internal  energy  in  the  gas  in  the  gun  tube,  joules. 

Kinetic  energy  of  the  gas  in  the  gun  tube,  joules. 

Kinetic  energy  of  the  liquid  in  the  gun  tube,  joules. 

Kinetic  energy  of  the  projectile,  joules. 

Kinetic  energy  of  the  piston,  joules. 

Conversion  constant  ■  10^  gm/s-cm-MPa. 

Liquid  enthalpy  in  the  reservoir,  joules/gm. 

Liquid  enthalpy  in  the  combustion  chamber,  joules/gm. 

Gas  enthalpy  in  the  combustion  chamber,  joules/gm. 

Liquid  enthalpy  in  the  gun  tube,  joules/gm. 

Gas  enthalpy  in  the  gun  tube,  joules/gm. 

Bulk  modulus,  MPa. 

Bulk  modulus  at  zero  pressure,  MPa. 

Derivative  of  the  bulk  modulus,  MPa. 


Gas  mass  in  the  combustion  chamber,  gm. 

Gas  mass  in  the  gun  tube,  gm. 

Liquid  mass  in  the  reservoir,  gm. 

Liquid  mass  in  the  combustion  chamber,  gm. 

Liquid  mass  in  the  gun  tube,  gm. 

Total  propellant  mass  in  the  system,  gm. 

Projectile  mass,  gm. 

Piston  mass,  gm. 

Mass  flux  into  the  combustion  chamber,  gm/s. 

Mass  flux  into  the  gun  tube,  gm/s. 

Rate  of  gas  formation  in  the  combustion  chamber,  gm/s 
Rate  of  gas  formation  in  the  gun  tube,  gm/s. 

Pressure  in  the  liquid  reservoir,  MPa. 

Pressure  in  the  combustion  chamber,  MPa. 

Average  pressure  in  the  gun  tube,  MPa. 

Pressure  at  the  gun  throat,  MPa. 

Pressure  at  the  base  of  the  projectile,  MPa. 
Projectile  resistance  pressure,  MPa. 
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Universal  gas  constant  =  8.318  joules/mole-K 


Projectile  travel,  cm. 

Piston  travel,  cm. 

Temperature  in  the  combustion  chamber,  K. 

Average  temperature  in  the  gun  tube,  K. 

■1 

Volume  of  the  liquid  reservoir,  cm  . 

•5 

Volume  of  the  combustion  chamber,  cm  . 

Volume  of  the  gun  tube  behind  the  projectile,  cm 

•3 

Liquid  volume  in  the  combustion  chamber,  cm  . 

o 

Gas  volume  in  the  combustion  chamber,  cm  . 

o 

Liquid  volume  in  the  gun  tube,  cm"^. 

■a 

Gas  volume  in  the  gun  tube,  cm  . 

Projectile  velocity,  cm/s. 
piston  velocity,  cm/s. 

Molecular  weight  of  the  gas,  gm. 

Porosity  of  the  combustion  chamber. 

Porosity  of  the  gun  tube. 

Ratio  of  specific  heats. 
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Liquid  density  in  the  reservoir,  gm/cm. 

Liquid  density  in  the  combustion  chamber,  gm/cm 
Liquid  density  in  the  gun  tube,  gm/cm. 

Gas  density  in  the  combustion  chamber,  gm/cm. 
Gas  density  in  the  gun  tube,  gm/cm. 


DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


No.  of 
Copies 


Organization 


2  Commander 

Defense  Technical  Info  Center 
ATTN;  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 


Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN;  SMCAR-ESP-L 
Rock  Island,  IL  61299 


1  - 


Director 

Defense  Advanced  Research 
Projects  Agency 
ATTN;  H.  Fair 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN;  AMSAV-E 
4300  Goodfellow  Blvd. 

St,  Louis,  MO  63120 


1  HQDA 

DAMA-ART-M 

Washington,  DC  20310 

1  Commander 

US  Army  Materiel  Command 
ATTN;  AMCDR-ST 
5001  Elsenhower  Avenue 
Alexandria,  VA  22333-0001 

3  Commander 

Armament  R&D  Center 
US  Army  AMCCOM 
ATTN;  SMCAR-TSS 

SMCAR-SCA,  M.  Devine 
SMCAR-SCA-T , 

R.  Yalamanchlll 
Dover,  NJ  07801 

5  Commander 

Armament  R&D  Center 

US  Army  AMCCOM 

ATTN;  SMCAR-LCA,  D.  Downs 

A.  Beardell 
SMCAR-LCE,  N.  Slagg 
SMCAR-LCS,  W.  Quine 
SMCAR-TDC,  D.  Gyorog 
Dover,  NJ  07801 

2  Director 

Benet  Weapons  Laboratory 
Armament  R&D  Center 
US  Army  AMCCOM 
ATTN;  SMCAR-LCB-TL 
P.  Votis 

Watervllet,  NY  12189 


1  Director 

US  Army  Air  Mobility  Rsch. 

and  Development  Lab. 

Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Communications 
Electronics  Command 
ATTN;  AMSEL-ED 
Fort  Monmouth,  NJ  07703 

1  Commander 

ERADCOM  Technical  Library 
ATTN;  STET-L 

Ft.  Monmouth,  NJ  07703-5301 

1  Commander 

US  Army  Harry  Diamond  Labs 
ATTN;  DELHD-TA-L 
2800  Powder  Mill  Rd 
Adelphi,  MD  20783 

1  Commander 

US  Army  Missile  Command 

ATTN;  AMSMI-R 

Redstone  Arsenal,  AL  35898 

1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-YDL 
Redstone  Arsenal,  AL  35898 


155 


DISTkIbUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


1  Commander 

US  Army  Belvoir  RiiD  Ctr 
ATTN:  STRBE-WC 

Tech  Library  (Vault)  B-315 
Fort  Belvoir,  VA  22060-5606 

1  -  Commander 

US  Army  Tank  Automotive  Cmd 
ATTN:  AMSTA-TSL 
Warren,  MI  48090 

Commander 

1  US  Army  Research  Office 

ATTN:  Tech  Library 
P.O.  Box  12211 
Research  Triangle  Park,NC 
27709-2211 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN :  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 

1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Banning,  GA  31905 

1  Commander 

US  Army  Development  and 
Employment  Agency 
ATTN:  MODE-TED-SAB 
Fort  Lewis,  WA  98433 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G33,  J.  East 
Dahlgren,  VA  22448 


2  Commander 

US  Naval  Surface  Weapons  Ctr. 
ATTN:  0.  Dengel 

K.  Thorsted 

Silver  Spring,  MD  20910 

1  Commander 

Naval  Weapons  Center 
China  Lake,  CA  93555 

1  Commander 

Naval  Ordnance  Station 
ATTN:  C.  Dale 
Code  5251 

Indian  Head,  MD  20640 

1  Superintendent 

Naval  Postgraduate  School 

Dept  of  Mechanical  Eng. 

ATTN:  Code  1424,  Library 
Monterey,  CA  93943 

2  AFWl/SVL 

Kirtland  AFB,  NW  87117 

1  Air  Force  Armament  Lab 

ATTN:  AFATL/DLODL 
Eglin,  AFB,  FL  32542-5000 

1  AFOSR/NA  (L.  Caveny) 

Bldg.  410 

Bolling  AFB,  DC  20332 

1  US  Bureau  of  Mines 

ATTN:  R.A.  Watson 
4800  Forbes  Street 
Pittsburgh,  PA  15213 

1  Director 

Los  Alamos  Scientific  Lab 
ATTN:  D.  Butler 
P.O.  Box  1663 
Los  Alamos,  NM  87545 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 


No.  of 
Copies 


1  Director  2 

Jet  Propulsion  Lab 

ATTN:  Tech  Lib 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 

2  Director  1 

National  Aeronautics  and 

Space  Administration 
ATTN;  MS-603,  Tech  Lib 

MS-86,  Dr.  Povinelli 
21000  Brookpark  Road  2 

Lewis  Research  Center 
Cleveland,  OH  44135 

1  Director 

National  Aeronautics  and 

Space  Administration  j 

Manned  Spacecraft  Center 
Houston,  TX  77058 

1  The  BDM  Corporation 

ATTN:  Dr.  T.P.  Goddard 
P.O.  Box  2019  1 

2600  Cearden  Rd,  North  Bldg 
Monterey,  CA  93940 

1  Calspan  Corporation 

ATTN;  Tech  Library 

P.O.  Box  400  2 

Buffalo,  NY  14225 

1  Food  &  Machinery  Corporation 

Northern  Ordnance  Division 
ATTN:  J.  Oberg 
Columbia  Heights  Post  Office 
Minneapolis,  MN  55421 

3  General  Electric  Ord.  Sys  Dpt  1 

ATTN:  J.  Mandzy,  OP4 3-220 

R.E.  Mayer 
H.  West 

100  Plastics  Avenue 

Pittsfield,  MA  01201-3698  , 


Organization 

General  Electric  Company 
Armanent  Systems  Department 
ATTN:  E.  Ashley 
M.  Bulman 

Burlington,  VT  05401 

Science  Applications,  Inc. 
ATTN:  R.  Edelman 
23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

TRW  Electronics  4  Defense 
ATTN:  D.R.  Ausherman 
M.P.  Bronstein 
One  Space  Park,  134/9048 
Redondo  Beach,  CA  90278 

Veritay  Technology,  Inc. 

ATTN:  E.  B.  Fisher 
4845  Millersport  Highway, 

P.O.  Box  305 

East  Amherst,  NY  14051-0305 
Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  Univ. 

Johns  Hopkins  Road 
Laurel,  Md  20707 

Director 

Chemical  Propulsion  Info 
Agency 

The  Johns  Hopkins  Univ. 

ATTN:  T.  Christian 
Tech  Lib 

Johns  Hopkins  Road 
Laurel,  MD  20707 

Pennsylvania  State  University 
Dept,  of  Hechnical  Eng 
ATTN:  K.  Kuo 
University  Park,  PA  16802 

University  of  Michigan 
Gas  Dynamics  Laboratory 
Aerospace  Engineering  Bldg. 
ATTN:  G.  M.  Faeth 
Ann  Arbor,  MI  48109-2140 


DISTRIBUTION  LIST 


I  No.  of 

Copies 

i 

2 


I 


Organization 

Princeton  Combustion  Rsch 
Laboratories,  Inc. 

ATTN:  N.A.  Messina 

M.  Summer field 
475  US  Highway  One  North 
Monmouth  Junction,  NJ  08852 

University  of  Mississippi 
Mechanical  Engineering 
Department 

ATTN:  C.R.  Wimberly 
University,  MS  38677 


1  University  of  Arkansas 
Department  of  Chemical 
Engineering 

227  Engineering  Building 
Fayetteville,  AR  72701 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 


Cdr,  USATECOM 

ATTN:  AMSTE-TO-F 

CDR,  CRDC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comment s/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ ^Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


CURRENT 

ADDRESS 


Name 

Organization 

Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


OLD 

ADDRESS 


Name 

Organization 

Address 

City,  State,  Zip 


(Remove  this  sheet  along  the  perforation,  fold  as  indicated,  staple  or  tape 
closed,  and  mail.) 


—  -  -  -  -  - - - FOLD  HERE  - 


Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground •  MD  21005-5066 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE.  taoo 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  WO  12062  WASHINGTON, DC 
postage  WILL  BE  MIO  BY  OERARTMENT  OF  THE  ARMY 


Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN;  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


FOLD  HERE 


